Phase transitions in one dimension: are they all driven by domain walls? 



Nikos Theodorakopoulos 1,2 
1 Theoretical and Physical Chemistry Institute, 

National Hellenic Research Foundation, 
Vasileos Constantinou 48, 11635 Athens, Greece 
(— 2 Fachbereich Physik, Universitat Konstanz, 
f—} 78457 Konstanz, Germany 



in 



(N 



O 



I 



^1" 

o 



-a 



(Dated: February 2, 2008) 



Two known distinct examples of one-dimensional systems which are known to exhibit a phase 
transition are critically examined: (A) a lattice model with harmonic nearest-neighbor elastic in- 
. teractions and an on-site Morse potential, and (B) the ferromagnetic, spin 1/2 Ising model with 

long-range pair interactions varying as the inverse square of the distance between pairs. In both 
cases it can be shown that the domain wall configurations become entropically stable at, or very 
near, the critical temperature. This might provide a "positive" criterion for the occurrence of a 
' phase transition in one-dimensional systems. 
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I. INTRODUCTION 



Phase transitions do not generally occur in one-dimensional systems. Their absence has been codified under a 
variety of conditions. Van Hove [1] showed that, for particle systems with pair interactions of sufficiently short range, 
the free energy does not have any singularities. Landau [2] argued that a bistable system with finite interface energy 
will achieve an entropic gain by breaking up into a macroscopic number of domains, making macroscopic phase 
coexistence at any positive temperature impossible. His argument explains the absence of a phase transition in a 
"-^J \ large class of one-dimensional systems, which includes the short-range 5 = 1/2 Ising model, the 4> 4 model, etc. 

How does theory cope with the exceptions from the broad rule? Proofs of the van Hove type are in essence proofs 
of analyticity. In principle they can be extended, refined and thus provide sharper criteria for the absence of a phase 
transition [3]. Nonetheless, it is not straightforward to make a positive statement about a system where a phase 
transition actually occurs. On the other hand, Landau's approach turns out to be more potent, because it can be 
t-H , "inverted" to describe the exact pathology which is necessary to produce macroscopic phase coexistence. Instead of 
setting out to prove that a finite interface energy will lead to a macroscopic number of domains because of the entropic 
gain invoved, it is possible to look at a single, possibly pathological interface, or domain wall (DW), and examine 
■ under which conditions the entropic gain might lead to its spontaneous formation. The question at hand is whether 
the conditions for spontaneous formation of a single interface are the same as those for a thermodynamic instability 
of the low-temperature phase, or, equivalently, for macroscopic phase coexistence. If the conditions are the same, i.e. 
' if (i) spontaneous DW formation occurs at some temperature T*, (ii) a macroscopic instability / phase transition is 
known by independent means to occur at a temperature T c , and (iii) T* = T c , then it is not unreasonable to argue 
1 that the transition is driven by the formation of the DW. 

In this paper I will review how this can be done in two representative cases: (A) models of one-dimensional ther- 
Q modynamic instabilities (interfacial wetting, DNA denaturation), (B) the 5=1/2 ferromagnetic Ising model with 
inverse-square interaction. Case (A) will be illustrated by the Peyrard-Bishop [4] Hamiltonian (continuum and dis- 
crete) - although this should be regarded as representative of a much broader class of on-site potentials with a repulsive 
core, a stable minimum and a flat top. A proper understanding of how DW formation governs one-dimensional phase 
transitions includes a discussion of symmetric on-site potentials, where the pathology is "accidentally" lifted and a 
phase transition does not occur. In a sense this represents a special case of (A) and will be treated accordingly. Case 
(B) has been recognized for a long time as a key representative of marginal pathology; in particular, Thoulcss [5] has 
used a Landau-type of argument to predict the occurrence of long-range order at low temperatures; his argument how- 
ever does not lead to an estimate of the critical temperature. It turns out that a "literal" application of the single- DW 
criterion (spontaneous DW formation), which has been successful in the case of thermodynamic instabilities, yields 
an estimate of the critical temperature which agrees with state-of-the-art Monte Carlo calculations [6] . This simple 
result is apparently new; it motivates the question stated in the title of the paper. In view of the pivotal role played 
by Serge Aubry in making our community aware of the relevance of nonlinear excitations (domain walls, solitons etc.) 
to statics and dynamics of phase transitions [7-9], I believe this is a most appropriate place to discuss it. 

The paper is organized as follows: Section II briefly restates the original Landau argument and a variant which 
will be useful here. The next two sections deal with cases (A) and (B), respectively. Conclusions are summarized and 
discussed in the final section. 
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II. A VARIANT OF THE LANDAU ARGUMENT 
A. Landau's original argument [2] 

Given a bistable system with N sites and an interface energy E separating phases A and B, the number of interfaces 
M is such as to minimize the total free energy 

F(M) = ME - TS(M) (1) 

where 



5 = In 



m 



'"I IT i 



_M\(N -M)\ 

is the entropy; the second expression in (2) involves use of Stirling's formula. Minimization of F with respect to M 
yields the most probable value of M 

M = Ne~ E/T . (3) 

At any T > the system breaks up into a macroscopic number of domains; in other words, phase coexistence cannot 
occur at a macroscopic scale. 



B. Comments and variations 



An obvious way to circumvent the prohibition (3) is an interface (domain wall (DW)) energy which grows loga- 
rithmically with the system size. If E = a In TV, the number of domains is macroscopic only for T > a. At T = a 
macroscopic phase coexistence can occur. Unfortunately, such a dependence of E on TV is usually accompanied by 
edge effects, i.e. the energy of the DW is no more independent of its position. In order to deal with such cases, it 
is more convenient to focus on a single DW, situated at position j, and examine its thermal properties. The DW 
partition function is given by 

N-l 

Zdw = ^ Edw{3)IT ■ (4) 

I will make explicit use of this in the context of case (B) below. For the moment, note that this formulation recovers 
the original Landau result if E(j) = EMj; in that case Zdw = N exp(—E/T), and the DW free energy is 

F DW = E-TlnN . (5) 

It follows directly from (5) that a DW will form spontaneously at a temperature T* = E/lnN, i.e., as long as 
E is finite, no phase coexistence at any nonzero temperature is possible in the thermodynamic limit. Moreover, if 
E = a In TV, the DW free energy will vanish - and phase coexistence can occur - at T* = a, just as predicted by (3). 

The vanishing of the DW free energy at a nonzero temperature provides an alternative, "positive" criterion for 
phase coexistence. I will use - and test - it in the remainder of this paper. 



III. CASE A: THERMODYNAMIC INSTABILITIES 



A. Background: Definitions 



Consider the class of onc-dimcnsional Hamiltonians 



N 



H = zZ 



n=0 



y + ^(jM+i-2Mr + V(y n ) 



(6) 



where y n , p n are, respectively, the displacement and momentum of the nth particle, R is a coupling constant and 
V(y) is any potential with a repulsive core, a stable minimum and a flat top. A convenient representative of this class 
is the Morse potential 



V(y) = (l-e-yf . 



(7) 



All quantities described above are dimensionlcss. The model defined by Eqs. (6)-(7) has been proposed in a variety 
of physical contexts, including interfacial wetting [10] and DNA denaturation [4]. 
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B. Background: Thermodynamics 

The thermodynamic properties of (6) are described by the classical partition function, whose nontrivial, configura- 
tional part Z can be expressed in terms of the spectrum of the symmetric transfer integral (TI) operator 

dy' K(y,y')My') = e- e " /T My) , (8) 

J — OO 

where 

K{y, y 1 ) = e -(v'-y) 2 /2RT e -[ V (y)+V(y')}/2T (9) 
In the thermodynamic limit N — > oo, 

Z = Y,e~ m " /T (10) 
is dominated by the smallest eigenvalue €q. The free energy per site is 

/ = ~lnZ«e . (11) 

"Analyticity" proofs examine the spectrum of the TI operator. For example, if the potential V(y) is unbounded 
from above as y —* ±oo, it can be immediately shown [11] that the operator in (8) is of the Hilbert-Schmidt type, 
and therefore the spectrum is nondegenerate at any T > 0, i.e. there are no singularities in the free energy. This 
scenario does not hold for the Morse potential - or any other potential of its class (cf. above) - which leads to a 
non-Hilbert-Schmidt kernel [12]. Finite-size scaling [13] suggests that a single bound state gradually merges into the 
continuum, i.e. 

eo(r)cx [T C {R)-Tf , (12) 

where T C (R) must, in general, be determined numerically. If the coupling constant is small, R -C 1, the harmonic 
coupling in (6) is strong. It is possible to make a continuum aproximation y n — > y(x). To leading order, this procedure 
in effect approximates the integral equation (8) by the Schroedinger-like equation 



RT 2 d 2 
2 dy< 



+ V(y) 



K{y) = tu<t> v {y) ■ (13) 



In the case of (13), the above scenario, i.e. the quadratic disappearance of the bound state into the continuum, can 
be verified analytically: Eq. 12 holds exactly, with 

t c (r) = Tr l (R) = 2 (|) 7 • ( 14 ) 

Other thermodynamic properties can also be calculated exactly; for example, it is possible to show [14] that the 
average displacement 

diverges linearly in the neighborhood of T c . 



C. The exact DW and its properties 



The equations of motion which correspond to the continuum limit of the Hamiltonian (6) admit exact static solutions 
of the type [14] 



y{x) = In 



1 + e 



±(a;— xo)/d 



(16) 



4 



where d = l/(2i?) 1 / 2 and xq is an arbitrary constant. These solutions have the property that they vanish as x — > ^oo 
and become linearly unbounded as x — > ±00. They can be regarded as interpolations from the bound to the unbound 
phase, i.e. they represent bona fide DWs. For any given transverse displacement of the boundary atom i/n+i = Y , 
yo = 0, the excess energy of the DW, compared to the minimum of the Morse potential, is equal to 2v, where 
v = Y/(2R) 1 / 2 is the number of unbound sites (not necessarily an integer, but this is not important for the present 
argument in the continuum limit). Their excess entropy, compared to the entropy of the bound chain, has been found 
[14] to be (R/2) 1 / 2 v. The DW excess free energy, 



Fdw — 



2-T 



1/2' 



(17) 



vanishes at T* = 2(2/R) 1 ' 2 , i.e. at the critical temperature (14). The above argument has been generalized [15] for 
DWs of the highly discrete (R ^> 1) version of the Hamiltonian (6); the numerical agreement between T* and T c is 
better than 1%. 



D. Symmetric on-site potentials 

It is well known that Schrocdinger-like equations of the type (13) with symmetric on-site potentials which are 
bounded from above will support bound states for any value of the control parameter R. The average value of the 
displacement < y > vanishes identically, since the ground state is of even parity. The system remains bound at all 
temperatures. 

The DW picture provides a straightforward explanation why such an on-site potential cannot generate a phase 
transition: in addition to the DW solutions (16) interpolating from the stable minimum to y — > +00, there is a 
symmetric pair of solutions interpolating from the stable minimum to y — > —00. Thus, although a particular DW 
solution can become entropically stable, an open system will always have an equal chance of creating the symmetric 
DW; the average displacement will vanish at all temperatures. 

IV. CASE B: THE FERROMAGNETIC S = 1/2 ISING MODEL WITH INVERSE SQUARE 

INTERACTIONS 

A. Background: Definitions, a glimpse of history & known estimates of T c 

The one-dimeensional, ferromagnetic S = 1/2 Ising Hamiltonian with long range interactions is defined by 

H= JyWn (18) 

where J > 0, <j„ = ±1 is known to have a phase transition if A < 2; the marginal case A = 2 is of particular interest; 
Thouless [5] predicted the occurrence of long-range order at low temperatures on the basis of a Landau-like argument. 
Anderson and Yuwal showed its relationship to the Kondo Hamiltonian and estimated K c = J/T c = 0.635 [16]. In 
spite of considerable theoretical progress during the years that followed [17], the exact value of the critical point 
remains unknown. A state-of-the-art Monte Carlo calculation by Luijten [6] gives an estimate K c = 0.6552(2). Other 
published estimates differ from each other by as much as 30% [6, 18]. 



B. The DW and its energy 

The spin configuration of (18) with minimal energy is cr„ = =Fl,Vn. A DW positioned at n = L is defined as a spin 
configuration with 

I will choose the upper sign; the excess energy Ejjw(L) of the DW (i.e. after subtracting the energy of the ground 
state, cf. above) can easily be calculated as 

E ™j L) = Sl(A - 1) - LTil(X) + £;v-l-i(A — 1) — (N — L)S w _ i _ 1 (A) - [Sjv_i(A - 1) - A^iv-i(A)] , (20) 
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FIG. 1: The excess energy of a DW positioned at the Lth site of a chain with N = 512, in units of J. The points represent 
exact values according to (24), the continuous line is the approximation of neglecting the contribution of f(L) (cf. text). 



where 

n , 

s «( A ) = E|x • ( 21 ) 

1=1 

Eq. (20) is valid for L < N/2; energies at larger values of L can be computed by using the symmetry Eow{L) = 
E DW (N~L). 

For A > 2 it is straightforward to establish that the DW energies are finite and essentially L— independent - except 
for a narrow range near the edges, which becomes vanishingly small in the thermodynamic limit. Eqs. (5) and (3) 
are both strictly applicable; either of them can be used to predict the absence of a phase transition. 

The marginal case A = 2 is more interesting. One needs 

" 1 

E »(2) = =C(2W(l + n) (22) 

i=i 

n 
1=1 

where 7 is the Catalan constant and ip, ip' the digamma and trigamma functions, respectively. It is then possible to 
rewrite the DW energy - now for the special case A = 2 as 



Edw{L) _ ^ 



N 



+ l + f(L) (24) 



2 J 

where the last term 

f{L) = tp(l + L) - Ini - 1 + Lip'(l + L) (25) 

is always small. This can be seen in Fig. 1, where I plot the exact energy and the approximation of entirely neglecting 
/; / can be shown to be entirely negligible for large L (of order (1 + L)~ 2 )\ it is largest at L = 1, but even then its 
value is t/>(2) + f(2) - 1 = 0.06772. 
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C. The DW's entropy: analytical approximation 

Neglecting the contribution of the /-term (cf. previous subsection), it is possible to write down an approximate 
partition function 

N-l r ^ IK 

(26) 



Dw ~ Mi(i-W 



where K — J/T. Three cases can be distinguished: 

1. high temperatures, K < 1/2: the sum in (26) can be approximated by an integral, where j/N is substituted by 
a continuous variable: 

2K 



Z lw = e-^«M>-™ J\ x {-^ ; ,27) 



in this case the DW free energy is given by 

i _ OK 

F DW = —\nN + 0(N°). (28) 

DW formation can proceed spontaneously; the change in free energy is of order — In N. 
2. K = 1/2: 

N-l 

7 * - p -(7+l) v - 

3 = 1 
N-l 



, 3(1 -j/N) 



T+D y fl + _LA 
f^\j N-jj 

N-l 

= 2e-^ +1 ) ~ 



3 = 1 " 

= 2e-^ +1) \nN + 0{l/N) . (29) 
DW formation can still proceed spontaneously, with a free energy change of order — In In N . 

3. K > 1/2: The leading term in the first half (j < N/2) of the discrete sum (26) - which is all that matters in the 
thermodynamic limit - can be obtained by neglecting the second factor in the denominator. Since the terms are 
symmetric around j = N/2, I estimate 

Z*dw « e~ 2 ^ K 2((2K ) . (30) 

Note that in this case Z* DW -and the DW free energy — h\Z* DW /K are of order TV . Spontaneous DW formation 
occurs according to whether or not Z* DW < 1. The critical temperature T* is obtained from the condition 
Z* DW = 1, i.e. 

^fP = l + T - (3D 

Numerical solution of (31) gives a value K* = 0.65136. Note that the value is quite close to the Anderson- Yuwal 
[16] estimate K c = 0.635. In fact, for K* not too far from 1/2, it is possible to use the limiting form of the £ 
function Q(2x) w l/(2a; — 1); in this case (31) becomes 

\n(K* - 1/2) , N 

- 2K« =1 + ^ ' (32) 

which is exactly the condition of [16]. 
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FIG. 2: The numerically calculated free energy of a DW in units of J. as a function of the coupling constant K. 
different form of iV-scaling for in the low and high K regime (cf. (28) and (30)). The dotted line marks the zero. 
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FIG. 3: The zeroes of the numerically calculated DW free energy plotted as a function of 1/lniV. The continuous line is a fit 
of the type (33) and yields K* = 0.6492(2). 



D. The DW entropy: numerical calculations 



It is possible to use the exact expression (24) of the DW energy and calculate the partition function (4) numerically. 
The resulting DW free energy as a function of K for a range of chain sizes is shown in Fig. 2. Note the different form 
of JV-scaling for in the low and high K regimes, cf. (28) and (30). 

The critical point is defined by the zero of the DW excess free energy (intersection of each curve with the dotted 
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line in Fig. 2.) A third-order polynomial fit of zeroes obtained for sizes 10 < TV < 8 x 10 7 , i.e. 

(cf. Fig. 3) yields the estimate K* = 0.6492(2) in the thermodynamic limit. 

V. CONCLUSIONS 

I have reviewed recent work on one-dimensional lattice models with asymmetric on-site potentials, which are known 
to exhibit thermodynamic instabilities, i.e. transitions from a bound to an unbound state, at a temperature T c . Such 
nonlinear models generically admit exact, unbounded solutions which can be regarded as domain walls. Examination 
of DW thermodynamics shows that a DW becomes entropically stable (i.e. the energy cost of producing it is balanced 
by the entropic gain) at a temperature T*. In the case treated analytically, T* = T C ; in the case treated numerically, 
equality holds within less than 1%. 

A similar procedure was successfully applied to the one-dimensional ferromagnetic S — 1/2 Ising model with inverse 
square interactions. In that case, the DW thermodynamics leads to a T* which also differs by less than 1% from the 
best available numerical estimate of the critical temperature. 

In summary, it appears that the vanishing of the DW free energy might be used as a "positive" criterion for the 
occurrence of a phase transition in one-dimensional systems, and at the same time provide an estimate of the critical 
temperature. In this sense, one might be tempted to conclude that phase transitions in one dimensional systems are 
indeed driven by DW formation. 
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